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The quantum mechanical two - particle problem is considered in hot dense nuclear matter under 
the influence of a strong electric field such as the field of the residual nucleus in heavy - ion reactions. 
A generalized Galitskii-Bethe-Salpeter equation is derived and solved which includes retardation and 
field effects. Compared with the in-medium properties in the zero-field case, bound states are turned 
into resonances and the scattering phase shifts are modified. Four effects are observed due to the 
applied field: (i) A suppression of the Pauli-blocking below nuclear matter densities, (ii) the onset of 
pairing occurs already at higher temperatures due to the field, (iii) a field dependent finite lifetime 
of deuterons and (iv) the imaginary part of the quasiparticle self-energy changes its sign for special 
values of density and temperatures indicating a phase instability. The latter effect may influence 
the fragmentation processes. The lifetime of deuterons in a strong Coulomb field is given explicitly. 
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I. INTRODUCTION 



The problem of Coulomb corrections in heavy - ion reactions raised a considerably interest during the last years. 
Especially, great attention is paid to the search of signals from the early stage of heavy-ion fragmentation. One of the 
most intensively studied experimental values is the two - particle correlation function. This observable gives hope to 
extract information about the process of fragmentation and the situation during the early stage of heavy-ion collisions. 
Coulomb effects from the residual nucleus are still not well investigated. Here we try to show that this effect may lead 
to new observable effects. Besides the strong interaction which acts exclusively at short distances of the fragments, 
the Coulomb field of the residual nucleus is a long range potential and will infiuence the evolution of the disintegrating 
nucleus for a longer time interval than the strong interaction. 

The most accepted way of describing heavy - ion reactions is to determine transport properties by solving adequate 
kinetic equations. One of the main ingredients of these equations are the transport cross sections. While a number 
of works have investigated the density dependence of the cross section and therefore the infiuence of the medium [1], 
modification to scattering properties due to the large electric fields have not been considered up to now. In this paper 
we estimate the infiuence of the Coulomb field of the residual nucleus the two-particle in-medium properties. The 
results are of interest for the transport theory of hot expanding nuclear matter in heavy ion collisions. 

A pioneering work in the area of field - dependent transport was done in solid state physics [2-7]. One of the 
problems is the derivation of correct transport equations including the field to arbitrary strengths [3,8-11,7,12]. The 
application of these high - field phenomena ranges from nonlinear conductivity [13-15] up to quantum coherence 
phenomena [16]. Among these studies there can be observed two different kinds of effects. The first class consists 
of kinetic effects resulting from the modification of the single particle distribution function in a strong field. As an 
example we refer to the stationary nonequilibrium function in strong electric fields [14-19]. A second class of effects 
are the quantum interference effects [13,20,21]. The latter one is considered in this paper. 

A very rough criterion for kinetic field effects is the ratio between energy gain of the charged particle during the 
mean free path A and the average energy lost during a collision. High field transport is characterized by an increase 
of the velocity between two successive collisions greater than the mean thermal velocity > Vt- The other class 

of field effects arise by quantum interference. These can be described by a typical time scale te = j^, which is 

determined for one particle properties by Se = ^( ^ sm^ Y^^. 

For two-particle properties [17] we will show that the coherence time is 



(^(£l + £i))i/3- 

V 8 V mi ^ mj 

If this time, which can be considered as a collisional broadening time, interferes with the time between two successive 
collisions one expects new effects of quantum nature [20,21]. The coherence time (1) can be seen from the two- 
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particle phase factor associated with times tl and t2. In the cm for the canonical momentum the phase reads with 
T = iti+t2)/2 and r = ti -^2 



V \\ 2n 2 \mi 1712 J J \n7i 7712 J 8 \77i1 7712 J 6 J J 

The term linear in E vanishes for equal e/m and the remaining quadratic term causes an intrinsic field effect. While 
the first term can be absorbed into redefinition of the energy scale the last term causes an intrinsic field effect with 
the parameter (1). This is immediately obvious by fouriertransforming to the frequency domain. 

Despite the general theoretical interest to describe two-particle properties in high electric fields we focus in this paper 
on the application to a system likely to be found in heavy ion reactions. To motivate the subsequent investigations we 
give a simple estimation of expected effects. During the expansion phase of highly excited nuclear matter after a heavy 
- ion collision the emitted charged particles will be subject to the Coulomb force of the remaining nucleonic system. 
For a nucleus of charge 100 one has a coherence time of about te ~ 100 fm/c at a distance of 10 fm. This becomes 
comparable with typical relaxation time scales in the expanding system. Considering the fact that this electric field 
will not change remarkably during the collision time we describe the effect of the radial Coulomb field in such a way 
that the spectator particle will get a short push by an effective constant electric field, where the time occurs only as 
a parameter. Instead of solving the complicated 3 body problem with correct Coulomb wave functions of the residual 
nucleus we propose to consider the transport during a short time period under the infiuence of a constant strong 
electric field. By this way we show how the two - particle properties will be modified owing to the high electric field 
impulses. This model treatment will be justified by the results discussed in chapter IV, in which we focus on three 
effects: (i) A suppression of the Pauli-blocking at densities no/40 — rio/lO, (ii) a pronounced onset of pairing states due 
to the field in comparison with the zero field case where it takes place at lower temperatures, (iii) the finite lifetime 
of deuterons due to the electric field and (iv) a modification of the density and temperature region of instability as 
indicated by the single particle self-energy. 

The outline of this paper is as follows. First we derive useful tools for formulating the theory gauge invariant. By 
the fundamental principle of gauge invariance we will present the in-medium T-matrix equation with the infiuence of 
field effects in chapter III. The dynamical behaviour is then discussed briefiy. In addition, we give a solution of these 
field- and medium- dependent Bethe-Salpeter equation for a separable potential of Yamaguchi-type. The interplay 
between field and medium effects is shown in chapter V, where we focus on the intrinsic field effect of quantum 
coherence. 



II. GREEN FUNCTION APPROACH 

We consider a system of charged nucleons under the infiuence of an applied constant electric field E. In vector 
potential gauge A(t) = — cEt, we have the Hamiltonian 

(ui - ^Ad))'-' 1 
H(t) = Y, + 2 (12|T^|l'2'>a+,a+,a2ai (3) 

1 ^ 1,2,1', 2' 

with a residual interaction V 

(12|T/|1'2'> = 6(h-t2)6(t[-t'2)6(h-t[)Vai(pi,P2,p'l,p'2)6p,+p,,p'^+p'^ 

which will be assumed separable. In order to describe high field effects and many-particle behavior correctly, it is 
necessary to use approximations for many particle effects which do not invalidate the resulting equations for arbitrary 
field strengths. It is known that the gradient expansion with respect to the times corresponds to a linearization 
in fields [11]. To get an unambiguous way of constructing approximations we have to formulate our theory gauge 
invariant. This can be done following a procedure known from field theory [22], which has been successfully applied 
to high field problems in [23]. First we present the idea for the one particle functions [12] and then for two- particle 
correlations. 

In order to describe correlations in highly nonequilibrium situations we define various correlation functions by 
different products of creation and annihilation operators in the Heisenberg picture 
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G>(l,2) = i<*i(ri,ti)*+(r2,t2)> 

G<(1,2) = < *^(r2,t2)*i(ri,ti) > . (4) 

Here < ... > denotes the average value with the unknown statistical noneqmUhnum operator p and =p signs the 
Fermi/Bose system respectively. 

A. Gauge invariant formulation of correlation functions 

As mentioned above, special care is required to formulate the theory gauge invariant. This will be achieved in the 
following way. First we change the space-time coordinates Xi = {ti,Vi) into Wigner coordinates x = xi — X2 and 
X = {xi + X'l)!'^. Then we construct gauge invariant correlation functions by a generalized Fourier transform [24] 

\ 

2 

I J d\ x^[k'' + ^A''(X + \x)] 

G{k,X)= j d'^x e -2 G{xX). (5) 



Here = ($, A) with $(R, T) and A(R, T) being the scalar and vector potentials, respectively. By gauge transfor- 
mations of j4'' and ^ one establishes immediately that G is gauge invariant. While it is possible to apply (5) to the 
general case, we restrict ourselves here to spatial homogeneous and steady electric fields. In this case (5) reduces to 
the simple rule [24] 

'4^ ^F[^M«^''+«rET], 



G{k,X) = J d'^x e^L""" +«"'^^JG(a;,X). 

Therefore, we must apply the following steps in order to formulate one particle correlation functions gauge-invariant: 

1. Fourier transform the difference- variable r into canonical momentum p 

2. Shift the canonical momentum p to kinematical momentum k according top = k — eET 

3. The gauge invariant functions G are then given by 

G(p, R, r, T) = G(k - eET, R, r, T) = G(k, R, r, T). (6) 

Next we generalize this method to two particle properties. Using the same Wigner coordinates as above, two-particle 
correlation functions can be written as 

Ga,{x,,x,,X,,X,) = + - + - ^)> (7) 

Within these coordinates it is possible to introduce the generalized Fourier transform (5) with respect to both difference 
coordinates 

Gabikl,k2,Xi,X2) = j d'^Xi j d'^X2 Gabixi, X2, Xi, X2) 

2 

X exp I J d\ (^(xi)^[k>^ + ^A''(Xi+\xi)] + (x2)^[k^ + ^A''(X2 + \x2)]). 

_ 1_ 

2 

(8) 

Consequently, the following rule is established to formulate two- particle gauge-invariant properties in the s-channel 
where ^1=^2, t[ = t'2, t = ti - t[ and T = (ti + t[)/2: 

1. Fourier transform the difference-coordinates rj = xj — Xj and r2 = X2 — X2 to the canonical momenta pi and 

P2 
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2. Shift the canonical momentum to kinematical momentum according to pi = ki — ejET 

3. The gauge invariant functions Gah are then given by 

Gah{pi,P2,'R-i,'R-2,T,T) = Gaj(ki - eiET, k2 - e-j'ET, Ri,R2,r, T) 

= G„i(ki,k2,Ri,R2,r,T). (9) 

In the subsequent considerations we will use yet another kind of coordinates which are introduced in the following 
way [25] 



< xiX2|G„,(tt )|x,x, >=J ^^^e 



dP dp dp' ^P(x,+x.-x^-x;,)/2 

l7r;i)3 (27r;i)3 (27r;i)3 
X eFP(^i-^2)+'FP'(^i-^2) < p|Gaj(P,R,tt')|p' > (10) 

where R = (xj + X2 + Xj + X2)/4 . Within these coordinates the gauge invariant form reads 

< p|G„i(P, R, T, t)\p' >=< p - £ii^ET|G„i(P - (ea + ei)ET, R, T, r)|p' + ^i^^ET > . (11) 
These rules are useful in formulating the two-particle properties gauge invariant and will be applied in chapter III. 

B. One-particle properties 

Next we give a relation between the different correlation functions of (4) which is not as obvious as in equilibrium case 
where we have the Kubo-Martin-Schwinger condition holds [26]. Assuming the conventional ansatz, i.e. replace the uj 
dependence of the distribution function by their quasi-particle value, the Wigner distribution function fw(p, R, T) = 
=FiG^(p, R,T = Q, T) can be related to the Green's function G^ by 

G<{pRuT) = TiaipRuT) fw(pRT). (12) 

This is quite good as long as the quasi-particle picture holds and no memory effects play a role. As we noted in chapter 
II A in the discussion following gauge invariance, the gradient expansion and the formulation of kinetic equations with 
high fields are basically connected with a careful formulation of retardation times. Therefore, the simple ansatz, called 
conventional Ansatz, will certainly fail. 

Another obscure discrepancy is the fact that with the conventional ansatz, one has some minor differences in the 
resulting collision integrals as compared with results from the density operator technique. With the conventional 
ansatz, one gets just one-half of all retardation times in the various time arguments [7,24]. This annoying discrepancy 
remained obscure until the recent work of Lipavsky, et al. [27], who established a modified ansatz. 

Before we use this ansatz we specify the spectral function including electric fields. For free particles corresponding 
to parabolic dispersions the gauge invariant spectral function under the infiuence of constant electric fields is well 
known [24,4] 

/ X „ 7 , 1 e2^2 
a^{p, w) = 2 / arcos— nuiT r r 



h \ 2m 24m 

b 

27r f p'^l2m-nuj\ , , 

Here Ai(x) is the Airy function [28] and we have defined ce = [Ti^ E"^ / I ^ . It is instructive to verify that (13) 
satisfies the frequency sum rule. In the subsequent considerations we will use this spectral function. The time variation 
of the field E{t) is assumed to be slow compared with Ti/eE- Then expression (13) can be used where E = E{T) 
depends on T as a parameter. 

This spectral function can be generalized to interacting systems. Within the quasiparticle picture the quasiparticle 
energy to = £p(R,T) is defined by the solution of the dispersion relation 

2 

ujh-- ReT.^{pujRT) = (14) 

2m 
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where is the retarded self-energy describing the many-particle influence. Then we flnd the field-free spectral func- 
tion as a sharp peak near the quasiparticle energies, which are now independent of to. Using Airy transformations [12] 
it is then possible to construct the held and medium dependent spectral function. In convenient time-representation 
one flnds [4] 



or in frequency representation 



asikTRT) = exp 



27r 

aEikwRT) = — Ai 



--[eURT)T+-—T' 
n \ 24m 



— {eu{RT) - ujh) 



(15) 



with k = p -|- eET according to (II A). This joint spectral function is a natural generalization of the quasiparticle 
picture to high-fleld situations and takes both effects into account: the gauge invariance and the many-particle 
influence. With the help of the joint spectral function (15) one obtains the ansatz valid for any applied electric held 
strength [12,17,27] in Hartree-Fock approximation 



GHkrRT) = Tie-(--^^1/. {k - £^,i?,T- 



(16) 



with the quasi particle energies e^- In order to get more physical insight into this ansatz one can transform to 
frequency representation [29] 

oo 2 2 

G<{kujRT) = Ti'i J dTcos^(^niUT-ekiR,T)T-^T''^ fwik-^,R,T-^). (17) 



Neglecting the retardation in fw one recovers the ordinary ansatz (12) with the spectral function (15) [see also 
(13)]. The generalized ansatz takes into account history by an additional memory. This ansatz is superior to that of 
Kadanoff-Baym in the case of high external flelds in several respects [21]: (i) it has the correct spectral properties, 
(ii) it is gauge invariant, (iii) it preserves causality, (iv) the quantum kinetic equations derived with Eq.(16) coincide 
with those obtained from the density matrix technique [7], and (v) it reproduces the Debye-Onsager relaxation effect 
[20]. 

III. T-MATRIX APPROXIMATION 

In order to describe short-ranged two-particle interactions it is necessary to introduce the standard approximation of 
the many-particle theory, the T-matrix approximation [30,25,26]. In the following we derive this approximation with 
complete time dependence and with the influence of arbitrarily high held strengths. This establishes a generalization 
of the known Bethe-Salpeter equation. Consequently, we flrst briefly repeat the general many particle derivation of 
the ladder approximation. This will be done carefully avoiding gradient expansions in time. In this way we obtain 
general expressions with complete time convolutions. Then we apply the derived transformation rules to formulate 
our equations gauge invariant. The Green's functions and the T matrix are given in gauge invariant forms. This will 
then be the stage at which approximations are discussed, since we still have reasonable control of held effects. 

Considering only binary-collision approximation the causal two-particle Green's functions can be written 



(121'2+) = G''(11')G'(22+) t (5„jG''(12+)G'(21') + i j dldVd2d2' 

[G''(11)G'(22) T'5aiG''(12)G'(21)] < 12|T„j|l'2' > G''(1'1')G'(2'2+) 



(18) 



Here the sum of ladder diagrams is deflned as a causal T-Matrix 

< 12|T„j|l'2' > = Kj(121'2') + i [ dld2d3d3'Vaiil233')G%31)G\3'2) < 12|T„j|l'2' > . (19) 



Since the interacting potential is assumed to be local in time, we can simplify the general equations. From the 
deflnition of the T-matrix (19) we get 
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With the abbreviation 



< 12|T|1'2' >=< xiX2ti\T\x[x'^t[ > 6{ti -t2)6{t[ - Q. 

< XiX2t\Q''^\xiX2i >= G''{xitXii)G^{x2tX2t) 



(20) 



we get the causal T-matrix 

< XiX2t\Tab\x'ix'2t' >= Vah(xi, X2, x[, x'2)6(t — t' ) 2 b{xi + X2 — x[ — x'^) 
+ i dxidx2dx3dx'^Vab{xi, X2, X3, x'^) < X2,x'3t\Qab\xiX2t >< XiX2t\Tab\x'ix'2t' > 



(21) 



Now we can apply the Lengreth-Wilkins rules [31,29] and find the relations in operator notation 



(22) 



By replacing the operator 1 — iVQ^ in the first equation with the help of the second one we derive the generalized 
optical theorem [30,32,33] 



(23) 



< xiX2t\T^\x\x'2t' > = j dxidx[dx2dx'2dtdt' < xiX2t\T^f,\xiX2t > 

> 

X < XiX2t\Q^fi\XiX2t' >< XiX2t' \T^\x'ix'2t' > . 

From the T-matrix equations (22) we get for the retarded one within the coordinates (10) 



< p\Tfi{PRtt')\p' >= Vai{p,p')6{t - t') + i 



d^p 



dtVai{p,p) < p\T^1{PRtt')\p' > 



G> (| - p, R, tt)G> (| + p, R, ti) -G<{^-p, R, tt)G<{^ + p, R, ti) 



(24) 



This system of equations was derived assuming small spatial micro-fiuctuations in comparison to macro-observables 
R. Thus, gradient expansions in space can be applied. It is important to remark that the equations up to now contain 
the exact time behavior for both microscopic times as well as center-of-mass times. Now we formulate these equations 
gauge invariant. 

We change the time variables according to t = t + f,T = (t + 1')/2, t = t — t' and apply the rules (11) to formulate 
the quantum mechanical two-particle problem gauge invariant. Then Eq. (24) becomes 



<p + 



^-^^ET\%^,{K, + Ki,R, T, r)\pi - '^±-^ET >= Vai{p,pi)6{T) 



d^p 



drVahiPyP) 



Gl{Ka-p+eaE- 



X G>(A'i +P + etE^^, R, -f, T + ^) - G<G< 



x{p- 



eh 



E{T + -)\fMiKa + Ki + (e„ + ei)E-,R,T+-,T + f)\p' 



(25) 



CaET and A'j 



eiiET have been introduced for compactness. 



The abbreviations Ka = 

Now we assume in agreement with our model of spatial and temporal homogeneous fields, that the macroscopic 
time T is much larger than the microscopic difference time r. This is justified if we consider the macroscopic time 
scales is in orders of the two- particle coherence time te = 100/m/c [see introduction and also the discussion after 
(29)]. Because r is smaller than r due to the integration limits, we see that T >> r and T >> r. In energy domain 
this means that we demand 
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For typical energies about 2bMeV this means that we restrict our considerations to macroscopic time T >> Afm/c. 
From this we expand r and r where they occur together with the macroscopic time T in Eq. (25) and obtain after 
introducing the ansatz (16) to zero order 

<P+ '-^^ET\%^,{Ka + Ki,T, t)\p' - '-^^ET >= Vai{p,p')6{T) 



+ if -(£^Vaiip,p) (1 - f^iKa-p,T) - fwiKi+P,T)) 

T 

X Jdx<p+ '-^^ET\ff^(Ka + Ki,x, T)\p' - '-^^ET > 



(26) 

Here we changed f = x — t. The energy arguments of the distribution functions are the quasi-particle energies. In 
appendix A and [17,12] a slightly more general result is found including retardation effects. Before we solve this 
equation in the next chapter we give the link to known results. In the field-free limit one recovers the standard result 
of the Bethe-Salpeter equation. To see this it is possible to write (26) in a more familiar way 



<P+'-^ET\ffjKa+K,,co,T)\p'-'-^ET>=Vai{p,p') + J J^Vai{p,p) 

.<p+ '-^ET\T-(K. + K„.,T)\p' - '-^ET > '^{^^ + - -) 

X (1 - (Ka -P,T)- fl ( A'i + p, T)) . (27) 
The function $ is given in terms of Airy-functions [28] 

<I>(z) = ^(^Gi(^) + iAi(^)y (28) 

Here the constant is A^; = Tive with the frequency ve introduced as 

h \ 8 nia mi, J 

As can be seen, the function $ reduces to l/(z + irj) in the field free case, where i>e tends to zero, such that the 
known Bethe-Salpeter equation occurs as the standard T-matrix equation [1,25,30,26,32]. Therefore we can interpret 
the inverse frequency te = Ve^ as the characteristic timescale describing the intra-collisional field effect. Similar time 
scales have been analyzed in detail in [34,35] for an electron-phonon system. The effect of finite collision duration 
is also presented in [36,6], where similar expressions for are discussed with respect to the band structure effects in a 
semiconductor. 

Without fields, the Bethe-Salpeter equation (27) with the factor l=p/=p/ = (1=F /)(1 T f) ~ f f includes, in 
principle, more than the Bloch-de Dominicis equation [37] with a respective factor (1 =p /)(1 =p /), by allowing for 
both intermediate particle - particle and hole - hole, excitations. In the zero-temperature limit for Fermions Eq. 
(27) corresponds to the Galitskii [38] equation, while the Bloch-de Dominicis equation corresponds to the Bethe 
- Goldstone [39] equation. Differences between the equations have been discussed in detail in the literature ( for 
references, see [30] ). 

IV. SOLUTION OF THE T-MATRIX EQUATION 

In consequence of the assumed separable potential of rank one 

Vab{p,p') = >^g{p)g{p') 
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we can represent the retarded T-matrix as 



(p|T„f (P, R, T, r)|p'> = Xgip' - ^i^ET)<,(p + ^i^ET)5(P - (ei + e2)ET, R, r, T). (30) 

Introducing this into eq (27) we can solve the equation for S by Fourier-transformation in difference time r. The 
solution of the T— matrix equation (26) reads 

/niT rp R ^ T)^\w\ - W-^Er)g(p + ^Er) 

(p|T„,(P,R,c.,T) |p)- i_AJ(P,R,E,c.,T) ^^^^ 

with 

J(P,R,E,c.,T) 

(32) 

The function $ is given by (28). In the previous equations the mass difference between the particle has been neglected. 
The more general expression is obvious. 

Equation (31) represents the main result of this paper. It gives the solution of the medium dependent Bethe-Salpeter 
or T-matrix equation in high fields including the medium effects via Pauli blocking factors. Similar investigations 
were performed in [7] within the framework of a random-layer model for semiconductors where the equations are less 
involved and the scattering T-matrix was used. 

Before we give numerical results in the next chapter, let us qualitatively discuss the physical content of (31). The 
field effects arise in two different ways. First, for different charged particles Ca 7^ ej, the applied field will decrease the 
T-matrix due to the coupling terms (ca — ej)ET and therefore will destroy any two particle correlation with increasing 
time. This is basically the fact that the two particles will be accelerated in different ways and lose coherence. In what 
follows we estimate the different time scales and show that the macroscopic time T is separated from the microscopic 
time r in the nuclear matter situation considered. The two particles become uncorrelated after a time T when the 
mean value of the momentum p < \ea — ei,\EcT corresponding to (31). Equating the latter expression we define the 
correlation time Tcorr = p/(\<^a — £h\Ec)- In a nuclear matter situation the effective radial Coulomb field of a nucleus 
with charge Z„ acting on a spectator is given by 

= f« % X 10« X V X fm. 

This means an elementary charged particle will feel a field of 100 MV/fm in a distance of 1 fm of the residual, 
charge-100 nucleous. The correlation time Tcorr is therefore Tcorr = 2'^(j^)^ for nuclear mater conditions. Now 
we return to the approximation which leads from eq (25) to (26). The macroscopic time T of the two- particle 
correlations are determined by the orders of Tcorr ■ The microscopic time r of eq. (25) is of the order of the time 
the particle needs to move the difference coordinate r. For this distance r the particle needs T = 3*^-^. Therefore 
the macroscopic time increases quadratically with the distance whereas the microscopic time is proportional to the 
distance from the residual nucleous. This means that our assumption of T >> r in deriving (26) from (25) is well 
justified. A generalization including retardation effects is given in the appendix A. 

For further numerical illustration of the content of equation (31) we will limit the discussion to equal charged 
particles like p-p correlations or assume for the triplet channel a sudden enough push of the electric field in order to 
establish stationary conditions - where only the intrinsic field effect is important. This second effect, which we focus 
on, is to be observed mainly for equal charged particles where the time dependence of the T matrix vanishes. The 
remaining coherence effect is condensed in the function $ in (28). In the field free case the function $ reduces to 
l/(z + i£) such that the known Bethe-Salpeter equation occurs [1]. This effect is a pure quantum mechanical one and 
is not possible to understand by classically. In the classical case one would expect that the correlation between two 
equally charged particles will not be changed by a constant electric field. Here we see that the correlations between 
two quantum objects can be altered by an applied electric field. This refiects the nonlocality of the scattering process. 
In semiconductor physics this effect is sometimes called intra-collisional field effect [24,34]. 
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V. NUMERICAL RESULTS 



For the purpose of exploratory calculation of the expected effects we choose a model for the interaction of Yamaguchi- 
type [40] 

Vab{p,p') = >^g{p)g{p') 

where 

^J(P) = 

The parameters A and /3 are chosen to reproduce the bound state energy of the deuteron and the scattering length 
for the triplet phase. 

In eq. (32) we performed the azimuthal integration and the remaining two dimensional integral was done numerically. 
It involves fast oscillating functions and is in general dependent on 3 angles, i.e. (p, E), (p' , E) and (A', E). In our case 
we consider the quasistationary condition, e.g. for equal charged particles, wherein the first two angular dependence 
vanish and only (A', E) remains. In the following we restrict to longitudinal fields with respect to the center of mass 
momentum K, which represents the maximal expected field effect. Due to the complicated field coupling a partial 
wave decomposition would be senseless here. Therefore we did not use angular averaged Pauli-blocking factors, but 
considered them completely. 



A. Scattering properties 

First we give the numerical solution of the T-matrix at positive energies. There the T-matrix becomes a complex 
valued function because J has a nonvanishing imaginary part. Similar to the low density limit, the on-shell scattering 
is characterized by phase shifts. Starting from the solution of the T-matrix in our model we consider the ratio 

tan h{p, n,T,E)= !™ ^ !I ^ , (33) 
xte < i > 

which carries certain information about the T -matrix and which has a simple interpretation as the momentum, 
density and temperature-dependent phase shift in the field free case [1]. In the considered field dependent case, this is 
a formal parameterization of the T matrix. The interpretation as phase shifts known from ordinary scattering theory 
are not obvious anymore. This is connected with the fact that one does not have plane waves as incoming particles. 
In high field case the incoming plane wave is modified to an Airy-function. The quantum mechanical scattering 
is therefore much more complicated. But the above definition can be understood as a quantitative measure of the 
medium and field dependent T -matrix and can be compared to the known phase shifts from field free solutions. 

For instance in figure 1 we plot the density dependent phase shifts for a temperature of 10 MeV versus the relative 
energy of two nucleons having a total momentum k = Q for the field- free case. All phase shifts were obtained with the 
numerical representation of the (5-function and the principal value, to which the expression (28) should collapse for 
vanishing fields. This highly oscillating function for small fields serves as a check for the numerics. The result agrees 
with the former ones in [41,1] as it should. For energies fc^/m < jj, the Pauli blocking factor (1 — — /^) becomes 
negative, so that this term in the Bethe-Salpeter equation acts like a repulsive force which gives rise to negative phase 
shifts at special densities [1]. This is to be seen for densities above rio/lO. 

From the behaviour of the phase shifts at A;^ = one concludes that in the low density limit the triplet state possesses 
a bound state following Levinson's theorem. As can be deduced from the figure this bound state vanishes between the 
densities no/20 and no/30 turning the phase shift from tt to 0. This behaviour will become more transparent when 
the bound states are considered explicitly in the next paragraph. 

For special poles Re(l — J{uj = 2//)) = the T-matrix approximation of the equation of state breaks down and the 
onset of ECS state can be observed [1,42]. In Fig. 4 the phase shifts are plotted for different temperatures at a fixed 
density. One sees clearly the enhancement of the resonance at fc^/m = 2// indicating the onset of superfiuity and the 
breakdown of the T— matrix approach. Further on we concentrate on temperatures above this critical value. 

In figure 2 we plot the same temperature and density dependent phase shifts but let an electric field act on the two 
particles. Here we introduce a scaling field-strength of 



xe-i X fm"^, (34) 
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where Za, Zi are the charges of the correlated particles. As can be seen we no longer have bound states as Levinson's 
theorem would suggest. This becomes more obvious in the figure 3 a where the low density phase shifts are plotted 
with and without fields. 

The second observation is that the Pauli-blocking effect becomes suppressed due to the field. Whereas in the field- 
free case the effective potential becomes repulsive at higher densities we see in figure 3 b that for densities around 
no/50 — no/20, where the field free shift already changes the sign, Pauli-blocking is less remarkable within electric 
fields. On account of the field additional phase space becomes available for scattering which is not accessible without 
fields essentially due to Pauli-blocking. This can be understood by the higly oscillating spectral function which leads 
to a collisional broadening resulting into a finite probability to scatter into higher energy states than the exact energy 
conservation would allow. This effect is not vanishing by averaging because of the highly asymmetric shape of the 
spectral function. Therefore it remains an additional probability to scatter beyond the Fermi energy due to the field. 

A third observation is that the onset of the singularity in the T-matrix due to the pairing is pronounced in a similar 
way as it is in the field-free situation, but occurs already at higher temperatures. As discussed in [1] the pole in the 
T-matrix, which is located at fc^/m = 2// in the field free and low density limit, gives rise to the onset of BCS states. 
The onset of this singularity is to be seen in the resonance around fc^/m = 2// which becomes more pronounced at 
low temperatures. In Fig. 4 we plot the phase shifts above and below the critical temperature for different densities. 
The sharp increase near fc^/m = 2// is demonstrated. The critical temperature is around 5 MeV. The same abrupt 
onset of transition to superfiuity is now observed with an electric field, but already at higher temperatures. Fig. 3 
c shows the phase shifts for a density no/10 and a temperature of T = 10 MeV where the field free shift changes 
sign and starts to show a resonance. In the field dependent case the resonance is very pronounced like it is for lower 
temperatures in the field free case. Therefore we believe that the pairing area in the temperature - density plane 
is enlarged due to field effects. This effect is very puzzling in the sense that from ordinary transport theory one 
expects an increase of the particle temperature and a destruction of coherence due to the field. Here we found quite 
the opposite, namely, for the two- particle problem the coherence phenomena leads to a higher pairing temperature, 
which would corresponds to lower kinetic temperature without coherence. This means the particles become ordered 
already at higher temperatures due to the field. This is understandable in the way that the coherence is built up. In 
other words, we obtain an enforced onset of coherence in a constant external field by the quantum mechanical nature 
of the scattering process. 

B. Bound states in high field 

Turning to the solution of the T-matrix with negative energies we can identify the poles as bound states in the 
system. In figure 5a we give the imaginary versus the real parts of the pole in dependence of the density at lOMeV. 
One sees that the bound state vanishes around no/20 — no /30. At these densities the imaginary part is zero translating 
into long living correlations, i.e. bound states. 

In the next figure 5b we give a plot of the imaginary versus the real parts of the poles at applied field strengths of 
20£'c corresponding to the scattering data in the last paragraph for the same varying values of density. We see that the 
poles acquire imaginary parts indicating the finite lifetime of the correlation. It is interesting to notice the nonlinear 
behaviour of the correlation lifetime with increasing density. First we see that with higher density the life-time 
increases slightly. At density regions where the Pauli-blocking effect will cause a change to a repulsive interaction, we 
see that the lifetime decreases with higher densities in the field-free case. The same nonlinear behaviour is observed 
in the next plot of real and imaginary parts vs. the applied electric field in Fig 6 for a low density limit. The energy 
of the correlation is reduced from the deuteron binding energy —2.22 MeV with increasing fields. At the same time 
the imaginary part rises sharply indicating that the bound state becomes a resonance. At a field strength of VbE^ 
the energy of the resonance turns back and crosses the real axes at a critical field strength of 20£'c- Correspondingly 
the imaginary part has a maximum indicating the shortest lifetime. The higher field behaviour is characterized by 
oscillations up to unstable behaviour, which is given by negative damping. 

The oscillating behaviour becomes more observable if one considers a density where the bound state is absent in 
the field free case. For a chosen density of no/4 in Fig 7 corresponding to the density in case d of Fig 1, we see 
that the pole moves in the complex plane crossing the real axes at higher critical fields than it is the case at lower 
densities. It is remarkable to observe that the damping in Fig. 7 is negative for small fields. This means that the 
system is unstable in situations where no bound states occur in field-free situations. The resonances are such that the 
correlations give rise to unstable behaviour in small fields, which may infiuence the fragmentation process. 
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VI. SUMMARY 



In the present paper we give have studied the problem of Coulomb signals from the early stage of heavy ion collisions. 
We have found that the influence of the Coulomb held of the residual nucleus on nucleonic correlations is important 
and can be estimated by a model of correlations which gets a Coulomb " push " in the early stage of evaporation or 
fragmentation. This modelled push of a constant held allows the formulation of the original three particle problem 
as a two-particle problem with an external held. By an explicit gauge invariant formulation of the nonequilibrium 
T-matrix we derive the Bethe-Salpeter equation with exact time evolution generalized to the influence on external 
electric fleld of arbitrary strength. 

Using a simple Yamaguchi-type of interaction we presented the solution of the Bethe-Salpeter equation with medium 
and fleld dependence. Field-free results reduce to results of earlier publications. An interesting interplay was found 
between medium and fleld effects. We observed that the two-particle bound states become resonances with corre- 
sponding fleld dependent lifetimes. A minimal lifetime occurs in dependence on fleld strengths. These bound states 
with flnite lifetimes are shifted towards positive values of energy which can be considered as the transition to scattering 
states. This indicates a scattering state as a two-particle correlation in the medium and fleld. 

Suppression of the phase shifts due to Pauli blocking is lowered by the applied fleld and can be understood in that 
the available phase space for scattering is enlarged by an applied external fleld. Furthermore, the onset of superfluidity 
as a resonance in the scattering phase shift is seen to become signiflcant already at higher temperatures in comparison 
to the fleld-free case. This is interpreted as the enforced onset of coherence due to the external fleld. 

By exploratory numerical solutions of the complex bound state energies we observe unstable behaviour at smaller 
fleld strengths and high densities which corresponds to the fleld-free situation where the bound states are already 
absent due to Mott transitions. 
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APPENDIX A: BETHE-SALPETER EQUATION WITH RETARDATION EFFECTS 

Here we like to present a less restrictive approximation to the exact equation (25), which was carried out in deriving 
eq. (26). It will lead to the inclusion of retardation effects. 

Instead of neglecting all difference times in (25) as compared to the macroscopic time T we now concentrate to the 
exact phase factor in (25). We introduce the ansatz (16) and the two particle propagator in (25) takes the following 
form 



exp I 



(Ka-p + e„E(^))2 ^ (Kb+p + eiE(^))2 



24 nia mi 



2ma 2mj 
X (^l-f^(K^-p-eaE(^+f),T+^ + f)-f^(Ki+p-eiE(^ + f),T+^ + f)) . (Al) 

If one now introduces (Al) into (25) one sees an complicated dependence on the difference time r as well as the 
integration time r. 

Besides the r dependence, which yields an oscillating function of Airy-type, there is an additional dependence on 
T+^^^^ in the arguments of the quasiparticle energies. This would yield additional oscillations besides the one already 
described. Therefore we restrict the attention to the mean oscillating function which was treated in chapter III. This 
means we should ensure that T + ^ T in the phase factor. This can be done in a saddle point like manner by 
replacing r — r in the remaining functions under the time integral. We obtain instead of (26) the slightly more 
general result 

<P+ '-^^ET\%^i{Ka + Ki,tT)\p' - '-^^ET >= Val{p,p')6{T) + i [ J^Vai{p,p) 



X I dx<p+ '-±^ET\f^^i(K^ + Ki + ^^i?(r), x,T- ^-)\p' - '-^ET > 
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X (l - f^(Ka -p- eaE^,T - ^) - /^(A'i +p- e,E^-,T 




(A2) 



In this result we account for retardation which leads to different dependence of the difference time r. The energy 
arguments of the distribution functions are the quasiparticle energies including off-shell effects by memory [29]. In the 
same manner as was done in deriving (27) we can express (A2) in a more familiar way using the Airy-transformation 



with p* = p -\- ^-^i-j^ET, p'* = p' — ^ ET and p* = p -\- ^ ET as abbreviations. Here the notation AT indicates 
the Airy-transformation 



where Ai(x) is the conventional Airy function [28] and the scaling frequency i>e is introduced in (29). As can be seen, 
the Airy transform will collapse to an identity transformation for zero fields, where i>e tends to zero, and we obtain 
the standard T-matrix equation provided we neglect the time retardation in the distribution functions. This result 
one would have been obtained without fields. In this way we derived a generalization of the Bethe-Salpeter equation 
including retardation effects in addition to the field dependence. 
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FIG. 1. The field-free triplet phase shifts for different densities vs. the energy difference of two correlated particles. For the 
low density limit the phase shift approaches the value tt indicating a bound state corresponding through the Levinson theorem. 
With increasing densities this bound state vanishes at a special critical density known as Mott transition. At higher densities 
the Pauli-blocking effects cause a sign change in the the phase shifts at small relative momenta equivalent to an effective 
repulsion. This leads to an onset of a resonance at energies around 2/t which can be interpreted as the onset of superfluidity. 

FIG. 2. The field-dependent triplet phase shifts for different densities vs. the energy difference of two correlated particles. 
The applied electric field strength are 20 x Ec- The same parameters are used as in Fig 1. The bound state at low densities 
turned into a resonance owing to the field effect and resulting in a value of the phase shift between and tt at small relative 
energies. 

FIG. 3. The field dependent triplet phase shifts (solid line) for 3 different densities together with the field free shifts (dotted 
line) vs the energy difference of two correlated particles. The applied electric field strength are 20 x Ec corresponding to Fig 
2. At small densities (top picture) the bound state is turned into a resonance resulting in a value of the phase shift between 
and TT. Further oscillations occur at higher energies due to field effects and the oscillating spectral function. At densities 
where the Pauli-blocking effects cause a repulsion in the field free case (middle picture) the field effects lead to a decreasing 
of Pauli-blocking. This is due to the enlargement of the available phase space by the external field. At high densities where 
the onset of the superfluid phase transition is starting as a resonance around energies of 2/t (bottom picture) the field effects 
cause an abrupt onset of the coherence. Without fields this behaviour is found at smaller temperatures, which means that the 
pairing occurs already at higher temperatures due to the field. 

FIG. 4. The field independent triplet phase shifts for 3 different temperatures chosen at two different density situations. 
The BCS resonance becomes enhanced with decreasing temperature. For nuclear matter density rio the critical temperature is 
about 5MeV. This illustrates the statement of the Fig 3 (bottom). 

FIG. 5. The imaginary versus real parts of the bound state energy of the deuteron state vs. density for a temperature of 10 
MeV for the field free case (above) and for a field of 20Ec (below). Whereas without fields (above) the bound state vanishes 
at a special density and the imaginary part is zero, at high fields the imaginary part is finite indicating a finite lifetime of the 
correlated state and shows a complicated density dependence of the resonance energy and damping. 
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FIG. 6. The imaginary (below) and the real part (above) of the correlation-state energy for vanishing density vs. the applied 
electric field. The vanishing damping or infinite life time of the deuteron becomes finite for higher field strengths and the bound 
state energy turns into a resonance moving in the complex plane. At special field strengths it crosses even the real axes and 
turns into an oscillating behaviour which is not plotted here. 

FIG. 7. The imaginary (below) and the real parts (above) of the correlation-state energy vs. the applied electric field for 
a density of rao/4 corresponding to the situation / of Fig 1. In contrast to Fig 6 a density is chosen where no deuterons are 
present without fields. One sees that the field dependence leads to a decreasing of the lifetime corresponding with an increasing 
of the imaginary part of the resonance energies at higher fields. It is pronounced that the small field behaviour shows a negative 
damping indicating an unstable situation. This is due to the fact that at densities of this high no deuterons are present anymore 
and the bulk nuclear matter becomes unstable against small perturbations in this density region. 
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